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Abstract 

The Gas Pixel Detector belongs to the very limited class of gas detectors opti- 
mized for the measurement of X-ray polarization in the emission of astrophysical 
sources. The choice of the mixture in which X-ray photons are absorbed and 
photoelectrons propagate, deeply affects both the energy range of the instrument 
and its performance in terms of gain, track dimension and ultimately, polarimet- 
ric sensitivity. Here we present the characterization of the Gas Pixel Detector 
with a 1 cm thick cell filled with dimethyl ether (DME) at 0.79 atm, selected 
among other mixtures for the very low diffusion coefficient. Almost completely 
polarized and monochromatic photons were produced at the calibration facility 
built at INAF/IASF-Rome exploiting Bragg diffraction at nearly 45 degrees. 
For the first time ever, we measured the modulation factor and the spectral 
capabilities of the instrument at energies as low as 2.0 keV, but also at 2.6 keV, 
3.7 keV, 4.0 keV, 5.2 keV and 7.8 keV. These measurements cover almost com- 
pletely the energy range of the instrument and allows to compare the sensitivity 
achieved with that of the standard mixture, composed of helium and DME. 

Keywords: X-rays, Gas Detectors, Polarimetry 
PACS: 29.40. Cs, 07.85.Fv, 95.55.Ka, 95.75.Hi 



1. Introduction 

Detectors able to image charged particle tracks in a gas have been developed 
over the last few years for different applications. One of the most promising 
is the possibility to resolve the path of photoelectrons emitted in the gas in 
consequence of a photoelectric absorption. The reconstruction of the initial 
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direction of photoelectron emission opens the way for measuring the state of 
polarization of the absorbed photons because the former is modulated with 
respect to the direction of the photon electric field with a cos 2 dependency. 
This makes the photoelectric effect a good analyzer of X-ray polarization, and 
a perfect one for absorption from spherically symmetric shells. 

Only a few gas detectors can resolve so finely the photoelectron tracks to 



accurately reconst ruct the initial direction of emission ([Bellazzini et all 12006 



Black et all 120071 ) . One of the most sensitive is the Gas Pixel Detector (GPP 



hereafter), developed by INFN-Pisa and INAF/I ASF-Rome (jCosta et all 12001 



Bellazzini et all 120071) and currentl y inserted in the focal plane o f several futur e 



satellite missions ( Bellazzini et all 12010 in presst ICosta et~aL , 2010 in press! ). 



The gas cell is 1 cm or 2 cm thick and a number of mixtures of helium, neon 
or argon and dimethyl ether (DME hereafter) at 1 atm or 2 atm have been 
used, the choice of the gas being of fundamental importance for the polarimet- 
ric performance of the detector. A hard limit to the lower energy threshold of 
the instrument is about twice the binding K-shell energy of the absorbing com- 
ponent because above this threshold the photoelectron track, modulated with 
polarization, prevails on the isotropic one of the Auger electron. Photoelec- 
tron range is determined by density, while the average atomic number fixes the 
mean free path for scatterings with atomic nuclei, which is the length scale on 
which polarimetric information is smeared. The diffusion coefficient influences 
the blurring of the photoelectron track during drift in the gas cell and even- 
tually the possibility to resolve the initial part of the photoelectron path and 
reconstruct correctly the direction of emission. We developed a Monte Carlo 
software to easily explore the behavior of the instrument to different mixtures 
and to subs equently test a subs et of the most interesting ones. 



Recently iMuleri et al.l (|2008l ) measured the modulation factor /i, namely the 
amplitude of the response of the instrument for completely polarized photons, 
for the GPD filled with helium 20% and DME 80% at 2.6 keV, 3.7 keV and 
5.2 keV. This data confirmed that measured values are basically consistent with 
what is expected on the basis of the Monte Carlo software and proved that X-ray 
polarimetry in Astrophysics with the GPD is feasible. In this paper we charac- 
terize the behavior of the GPD with a different gas, i.e. pure DME at 0.79 atm. 
In particular, we describe the configuration of the GPD and the calibration 
sources we used in Section [2j while spectral capabilities of the instrument are 
discussed in Section [3] The measurement of the modulation factor between 
2.0 keV and 7.8 keV is reported in SectionHJ together with the comparison with 
Monte Carlo results and what was reported previously on the He 20% and DME 
80% mixture. Note that this is the first time that the modulation factor of a 
gas polarimcter is presented at energies as low as 2.0 keV. 

2. Set-up 

2.1. Detector configuration 

The Gas Pixel Detector is composed of a sealed 1 cm thick gas c ell en 



closed by a 50 /J,m beryllium window, a GEM (Gas Electron Multiplier, ISauli 
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19971) which collects and amplifies primary electrons produce d by photoelectrons 



in the gas cell, and a finely subdivided pixelized detector ([Costa et all 12001 



Bellazzini et all 120071 ). The last component, based on a VLSI ASIC realized 
in 0.18 iim CMOS tec hnology, is the actual breakthrough of the instrument 
( Bell azzini et al ., 200G), which otherwise is fundamentally an array of standard 



yet exceptionally small independent proportional counters. The top metal layer 
of the CMOS is fully pixellated to collect the charge produced in the common 
gas volume and allows to obtain a true 2D image of the photoelectron track even 
at low energy, thanks to the small (50 fim) pixel size. The acquisition is self- 
triggered and only a small window of about a thousand of pixels enclosing the 
track is actually read-out in place of the whole matrix. The chip is 15 x 15 mm 2 
and comprises 105,600 pixels arranged in a hexagonal pattern. 

The cell is sealed but can be refilled to test different gases and typically 
mixtures of helium, neon or argon and DME are used. DME is used to reduce 
diffusion and also as a quencher, but it acts as the actual absorber in the case 
of helium mixtures. The first application of the instrument in Astrophysics is 
expected in the 2-10 keV energy rang e and within this i nterval the standard 
mixture is helium 20% and DME 80% (|Muleri et all 120081) . This was preferred 



to mixtures of neon because of the longer photoelectron path and lower diffusion 
for equivalent efficiency, which assure a higher polarimetric sensitivity at low 
energy where the largest part of photons are concentrated. In this paper we 
push the use of low-diffusion mixtures to the extreme, exploring the use of a 
pure DME gas at 0.79 atm (0.8 bar). Since helium is basically transparent to 
X-rays in the 2-10 keV energy range, we expect a sensitivity comparable to the 
standard mixture, with a possible enhancement because of the lower diffusion. 

An improvement with respect to previous versions of the GPD is the use of 
a laser-etched GEM made of liquid crystal p olymer which shows a better tem- 



a laser-etched (jr^lvi made or liquid crystal p olymer winch snows a better tem- 
poral gain stability ( Tamagawa et all 120091) . A drawback is that the smallest 



pitch available was only 80 /im (instead of 50 ^m of previous detectors), and 
this has proved to be insufficient to avoid the emergence of systematic effects 
due to undersampling of short tracks, discussed and removed, as explained in 
Section B~T1 GEMs with smaller pitch are now in production and will be used for 
the next GPD prototype. Moreover the thickness is 100 /im instead of 50 /im. 
The characteristics of the GPD used are summarized in Table [TJ 

2.2. Calibration facility 

The GPD was characterized at the X-ray facility built at INAF/IASF-Rome. 
Although its detailed description is beyond the scope of this paper, in the fol- 
lowing we briefly present what is relevant to measurements presented below. 

P olarized and mono chromatic X-rays are produced by Bragg diffraction at 
45° (|Evans et all Il977l) . Incident radiation on a crystal can be decomposed 



in two components, polarized parallel (7r-component) and perpendicularly (cr- 
component) to the diffraction plane. The latter is more effectively diffracted 
because the ratio k between the integrated reflectivity of the tt and a components 
is always smaller than 1. Hence diffracted radiation is (partially) polarized and 
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Area: 


15x15 mm 2 


Active area fill fraction: 


92% 


Window: 


50 /im, beryllium 


Mixture: 


DME 100%, 0.79 atm (0.8 bar) 


Cell thickness: 


1 cm 


GEM material: 


copper-coated liquid crystal polymer 


GEM pitch: 


80 fim 


GEM holes diameters: 


48 fiva 


GEM thickness: 


100 (jxa 


GEM voltages: 


V d „ /t =3200 V, V top =lU5 V, V 6o « om =500 V 


Gain: 


500 


Pixels: 


300x352, hexagonal pattern 


Pixel noise: 


50 electrons ENC 


Full-scale linear range: 


30000 electrons 


Peaking time: 


3-10 /is, externally adjustable 


Trigger mode: 


internal, external or self-trigger 


Self-trigger threshold: 


2000 electrons 


Pixel trigger mask: 


individual 



Table 1: Main characteristics of the GPD prototype studied in this paper. 



the degree of polarization V is: 



V = 



1 - k 



1 



(1) 



If the incident angle 9 is 45°, k — and consequently V = 1. For intermediate 
values, k can be calculated and the value as a function of 9 is reported in 
Figure 1(a) for graphite crystals. The large dependence of k on the incident 
angle requires the value of 9 to be tightly constrained to prevent the dilution 



of the average degree of polarization (cf. Figure 1(b) ). The angular constraint 
also selects the energy of diffracted radiation, related to 9 by Bragg's Law: 



E(9) = 



ihc 



2d sin 9' 



(2) 



where h and c are respectively Planck's constant and the speed of light, d the 
crystal lattice spacing and n an integer which specifies the diffraction order. 

We already presented a prototype source based on Bragg diffraction, which 
exploits lead-glass capillary plates to constrain to 45° the incident and diffraction 
angles and small (2 W) X-ray tubes to produce the radiation to be diffracted 
(jMuleri et all 20071 ). This source was used to generate polarized photons at 
2.6 keV, 3.7 keV and 5.2 keV and calibra te at these energies the GPD filled with 
a Hc-DME mixture (jMuleri et all 120081 ). An aluminum crystal and an X-ray 



tube with anode made of calcium were exploited to produce 3.7 keV polarized 
photons, while 2.6 keV and 5.2 keV were obtained by first and second order 
diffraction on graphite of copper X-ray tube radiation. The former configuration 
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Figure 1: (a) Dependence of k with incident angle and energy, relate d by Bragg's law, for 
graphite crystals. The value of k was calculated by iHenke et alj ||1993T ) (b) Expected degree 
of polarization, derived from Equation JTJ- 



is particularly effective in terms of higher flux and control of the output state of 
polarization because Ka fluorescence emission of calcium is well in accordance 
with Bragg energy at 45° for aluminum. Then (almost) all incident photons have 
a well-defined energy, that of Ka line of calcium, and are diffracted exactly at 
the Bragg angle given by Equation @ and the degree of polarization is precisely 
calculated with Equation ((T|). A trade-off between flux (low collimation) and 
high polarization (high collimation) was instead necessary for diffraction on 
graphite because X-ray tubes with anodes in accordance with Bragg energy are 
not available in this case and continuum bremsstrahlung emission is to be used. 

After the construction of this first prototype, we built a more powerful source 
based on the same concept but with some differences. The most important 
is that the crystal is mounted on a manual stage which allows two axes tilt 
regulation in the range ±3° to achieve the best alignment to the Bragg condition. 
More crystals are available to produce radiation at different energies and, thanks 
to more powerful X-ray sources (50 W) , a tight collimation can be retained even 
for diffraction of continuum radiation, which is less effective than the use of 
line in accordance to Bragg condition. A mechanical assembly combined with 
motorized and manual stages, which complete what we call X-ray facility, allows 
the detector to be moved, rotated and inclined with respect to the beam (see 
Figure^). 

Within the work presented in this paper, we used three different crystals, 
PET, graphite (grade D, 1.2° of mosaic spread) and aluminum. In the first 
two cases, we used the improved version of the Bragg source and exploited the 
diffraction of continuum photons produced by an X-ray tube with anode made 
of titanium, which has the thinnest window among our medium power tubes. 
We constrained both incident and diffraction angles on the crystal with two 
capillary plates which provide a collimation 1/40 (semiaperture 1.4°) and 1/100 
(0.6°) respectively. The PET crystal was used to generate 2.0 keV and 4.0 keV 
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Figure 2: Setup of the measurements. X-rays, generated by a 50 W tube manufactured by 
Oxford Instruments, are emitted horizontally and diffracted downwards by a crystal oriented at 
45°. Two collimators constrain the divergence of the incident and of the diffracted photons and 
a diaphragm, placed after the second collimator, limits the size of the beam. The detector is 
mounted on a platform which is allowed to rotate, move and incline with respect to the incident 
beam with high precision, < f (im for movements and < 1 arcsec for rotation/inclination. 
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Energy Crystal X-ray tube Collimation Diaphragm 

(keV) mm 



2.0 PET, 1st order 
4.0 PET, 2nd oder 



Ti, 50 W In: 4, Out: ^ 2 >0.99 



2.6 GrD, 1st order 
5.2 GrD, 2st order 
7.8 GrD, 3st order 



5.2 GrD, 2st order Ti, 50 W In: 40' 0ut: 100 02 >0.99 



3.7 Al, 1st order Ca, 2 W Out: ± 2 0.9938 

Table 2: Setup for each crystal. "GrD" stands for the mosaic graphite grade D crystal. 



polarized radiat ion, corresponding to the first and second order of diffraction 
{2d = 8.742 A, iHenke et alJ , Il993h . while for graphite w e exploited the firs t 
three orders at 2.6 keV, 5.2 keV and 7.8 keV (2d = 6.708 A. lHenke et"aT1 . ll993f ). 
At 3.7 keV we used the prototype source because of the favorable accordance 
between the aluminum crystal and the Ka line of the calcium tube. In this 
case the geometry was already constrained by the use of nearly monochromatic 
incident photons and hence only a single collimator (with collimation 1/40) 
limited the diffraction angle. Helium flowing in the Bragg diffractometer was 
required with PET crystal to avoid severe air absorption of 2.0 keV photons. 
A diaphragm 2 mm in diameter was used to illuminate only the small central 
region of the detector. The setup for each crystal is summarized in Table [5J 

Spectra diffracted by crystals were acquired by means of a Si-PIN Amptek 
XR100CR spectrometer with a energy resolution ~200 eV FWHM at 5.9 keV. 
Lines were fitted with a gaussian profile and the measured line energy was used 
to estimate the degree of polarization V by the calculation of the diffraction 
angle and of the interpolated value of k. Thanks to the tight collimation of 
the two capillary plates on the incident and diffraction directions, the degree 
of polarization is always above 99%. The lack of any significant dilution of 
polarization is also confirmed by the narrowness of diffracted lines, which have 
a FWHM (-200 eV, 215 eV at 7.8 keV) consistent with that of fluorescence 
lines at comparable energies. 

As an example, we report in Figure [3] the spectrum for diffraction on the 
graphite crystal of radiation generated by the titanium X-ray tube. The first 
three orders are prominent at 2.6 keV, 5.2 keV and 7.8 keV but from residuals 
is also evident a line at about 3.5 keV, which is the escape peak from silicon 
of the 5.2 keV photons, and a contribution at 4.5 keV. The latter is caused by 
Ka titanium fluorescence photons scattered on the crystal or on its holder and 
reaching the detector after passing through the two collimators. By the way, the 
two collimators are completely opaque at this energy and then only scatterings 
at nearly 90° are possible. This implies that even scattered photons are highly 
polarized perpendicularly to the scattering plane. Even if the GPD is not able 
to resolve this line and that at 5.2 keV, the flux of the former is < 1/300 than 
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20 




Energy (keV) 

Figure 3: Spectrum obtained by diffraction on the graphite crystal of continuum photons pro- 
duced with the titanium X-ray tube and acquired with the Amptek XR100CR spectrometer. 
The first three orders are visible at 2.6 keV, 5.2 keV and 7.8 keV and gaussian fits of each line 
are also reported as red solid lines. The high voltage of the X-ray tube was set to 10 kV to 
avoid the presence of photons at energies above the third order. The line at about 3.5 keV is 
the escape peak of 5.2 keV photons and that at 4.5 keV is the scattered titanium Kct. 

that at 5.2 keV and therefore it will be neglected in the following. 

3. Spectral capabilities 

In Figure 0] we report the spectrum acquired with the GPD for monochro- 
matic (and polarized) radiation obtained by Bragg diffraction on PET (2.0 keV 
and 4.0 keV) and graphite crystals (2.6 keV, 5.2 keV and 7.8 keV). Energy of 
photons is derived by the sum of the charge content of hit pixels. The lines are 
well resolved and this allows us to present in the next Section the modulation 
factor for each energy without any significant crosstalk. Lines are fitted with 
a Gaussian profile and the results are reported in Table [3] Data is filtered (i) 
to select spatially the incident beam and (ii) remove the tracks with two or 
more clusters of hit pixels, that is we postpone those events (<5%) which are 
composed of non-contiguous groups of pixels to subsequent refined analyses. 

The relation between the energy and the pulse height, reported in Figure[5J is 
fitted with a line y = mx + q, where q = (-239.7±2.3) ADC and m = (1055.8 ± 
0.8) ADC/keV. The deviations from linearity (5% at 7.8 keV) can be explained 
by the insufficient control of high voltage values and/or the inadequate monitor 
of the environment conditions. The HV power supply used (CAEN N470) has 
a accuracy ±1 V and a long term stability ±2 V, which imply variations on the 
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Energy 




FWHM 


5E/E 


(keV) 


(ADC) 


(ADC) 


(ADC) 


2.0 


1867 ± 1 


689 ± 3 


0.369 ± 0.002 


2.6 


2480 ± 1 


797 ± 3 


0.321 ± 0.001 


3.7 


3783 ± 2 


1071 ± 4 


0.283 ± 0.001 


4.0 


4013 ± 3 


1104 ± 8 


0.275 ± 0.002 


5.2 


5099 ± 3 


1248 ± 8 


0.245 ± 0.002 


7.8 


7645 ± 15 


1793 ± 23 


0.235 ± 0.003 



Table 3: Results of the fit to the lines at different energies acquired with the GPD. 
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Figure 5: Energy calibration of the GPD. 



gain of the order of ±4%. The gain G can be normalized by temperature T and 
pressure p, expressed in Kelvin and torr respectively, with the function: 

r^corr i 

T (3) 



f^meas 

^ exp 



Cx u 



Jp/T 2.533 y 

where G meas and G corr are the measured and corrected values and C is a con- 



stant, equal to 19.1 Torr/K for the RIKEN-80T-LCP GEM ([Tamagawa et al 



2009). Changes of temperature cause linear variations on the gain, Ag cx — AT, 



and fluctuations of 2 K, which are plausible in our case, results in gain instabil- 
ities of 5%. 

Within Poisson or Poisson-like statistics of ionization and charge amplifica- 
tion (Knoll, l2000h . the dependence of the energy resolution 8E/E = FWHM/£ p , 



CO k 



as a function of energy is expected to be 6E/E cx l/y/~E for a gas detector. The 
relation in the case of the GPD is reported in Figure [51 although the energy res- 
olution decreases with energy, the sole dependence with \j\fE is unsatisfactory. 
Therefore we used a function: 

+ (*a E)\ (4) 




where k x = (0.520±0.001) VkeV and fc 2 = (0.0188±0.0005) keV 1 (dashed line 
in Figure [6]). The fit function comprises two contributions summed quadrati- 
cally, the first being the statistical Poisson noise in the ionization and ampli- 
fication of the charges. We used the second term to model in a simple way 
an additional contribution which seems to weigh more at higher energy. A 
significant role within it could be played by incomplete charge collection due 
to photoelectron hits on the gas cell surfaces, which depends on the direction 
of emission but also on energy through the range (1 mm at 7.8 keV). We po- 
sitioned the beam in the center of the detector, but loss of charge can still 
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Figure 6: Energy resolution as a function of energy. The fit function (dashed line) takes 
into account the statistical noise and a linear contribution. These two terms are plotted with 
dotted lines. 



occur if photons are absorbed near the window or the GEM. More energetic 
and hence longer tracks are also spread on larger regions of the GEM and of 
the ASIC and then nonuniformities could give a contribution. While the for- 
mer have be en reported at the level of several percent on spatial scales of a few 
millimeters ( Tamaeawa et alll2009t l. we are confident that ASIC contribution 



is small if not negligible. Pedestals are read-out and substracted immediately 
after the event (delay o f a few /xs) and pixel noise is very low, 50 electrons ENC 
( Bcll azzini et al ., 200G). Tracks produced by photons at 2.0 keV hit on average 
40 pixels (see below) and then the total contribution of electronic noise is of the 
order of \/2 a/40 50 « 450 electrons, where the factor y/2 is for pedestals sub- 
straction. Conversely, 2.0 keV photons produce 71 primary electrons, assuming 
that t he average energy lo ss for the creation of a ion-electron pair in DME is 
28 eV (jPanskv et allll997l) . and Poisson fluctua tions are y/J 71, where / is the 
Fano factor -0.3 for DME (jPanskv et al.1 . Il997h . After the GEM amplification 
of a factor 5 00 (see Table [T]) , statistical fluctuations result in 4800 electrons 
(iKnolll . l2000h . i.e. a factor 10 larger than electronic noise. Actually, we expect 
that the main source of the possible electronic contribution to energy resolution 
is the spatial nonuniformities in pixels gain, which can't be calibrated because 
sufficient precise test capacitances does not allow to be build with the 0.18 mm 
VLSI technology. 

Spectral performance achieved is overall quite good. The energy resolution 
at 5.9 keV is less than 24% and then very close to the requirement of 20% 
for space missions that foresee a GPD on-board, which, however, assume that 
energy is measured with a dedicated spectroscopic channel taking the signal 
from the upper electrode of the GEM. 

Interestingly enough, the average number of hit pixels (p) is a very well 
defined function of energy (see Figure [7]). The dependency can be modelled 
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Figure 7: Average number of hit pixels per track p as a function of energy. Errors are evaluated 
as the ratio between the root mean square of the p distribution and the square root of counts. 



with a function: 

E \ k 



(P)=k a + k b , (5) 



1 



with k a = (31.15±0.13), k b = (3.58±0.05) and k c = (1.664±0.008). Note that 
the index of the power law is broadly consistent with the energy dependence 
of electron range i n a gas, rep orted to be R = 0.71 -EmcV g/cm 2 under a few 
hundred of keV by|Sauljl (1977J). 



In principle the Equation ([5]) could be exploited to derive an estimate of 
the photon energy, but unfortunately the correlation is rather weak because 
the number of hit pixels has also a strong dependence on the photoelectron 
direction of propagation. Photoelectrons emitted orthogonally to the collection 
plane produce smaller tracks than those which instead propagate mostly parallel. 
The errors associated to the energy determination by means of the Equation ([5]) 
could be derived precisely only by the knowledge of the probability function of 
p, but with an approximate gaussian fit the energy resolution is rather poor 
(> 60%, about constant with energy). 



4. Modulation factor 

In general, the response of a polarimeter is modulated as a result of the 
absorption of polarized photons: the degree and the angle of polarization are 
derived from the amplitude and the phase of this modulation, respectively. The 
modulation curve M(4>) of a photoelectric polarimeter is the number of pho- 
toelectrons emitted per each azimuthal angle 4>. The direction of emission is 
most probably aligned with the electric field of the absorbed photon, th e depen 



dency being expressed by the differential cross section of the interaction (jHeitler 



19541) : 

ph 



daH sin 2 9 cos 2 cf> 

« ^, (6) 



dn (i + /3< 
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where /3 is the photoelectron velocity in units of c. 9 and <j> are the latitudinal and 
azimuthal emission angles, measured with respect to the direction of incidence 
and of polarization respectively. The formula reported above is valid only in 
the case of spherically symmetric shells, otherw i se cor rections which reduce the 



response to polarization are required (jGhoshl . 119831 ). As far as the GPD is 



concerned, K-shell is largely the most probable involved in photoabsorption and 
then the Equation ([6]) is sufficiently accurate. 

For completely polarized photons, the number of photoclcctrons emitted in 
a certain direction is modulated as a cos 2 function, while for partially polarized 
radiation the amplitude is reduced linearly. Real instruments do not respond 
perfectly to polarization, in the case of the GPD because of scatterings with 
nuclei, gas diffusion and finite size of ASIC pixels. Hence even for 100% polarized 
photons the amplitude of the response, called modulation factor /x, is never 
complete but it is calculated as: 

^ M max + M III i II 

where M max and M. m in are the maximum and the minimum of the modulation, 
which practically is the histogram of photoelectron angles of emission, in the 
case of completely polarized photons. Since A4((f)) is fitted with a function 
M(4>) = A + B cos 2 ((f>- O ), then: 

The modulation factor, and the efficiency e, are the primary parameters used 
to derive the sensitivity of a polarimeter. The modulation measured, propor- 
tional to //, is to be compared with that naturally arising from statistical Poisson 
fluctuations of the number of photoelectrons emitted per angular bin, which are 
inversely proportional to y/e if the background is negligible as for an experiment 
at t he focus of an X-ra y telescope. The product \iyfi is called quality factor 
(see iMuleri et all (|2008l ) for a more extended discussion on the sensitivity of 



a polarimeter). Within astrophysical application, the value of n is especially 
important at low energy, where the largest part of photons are concentrated. 
Typical source spectra rapidly decrease with energy, as the efficiency of the in- 
strument and the area of grazing incidence optics. For a Crab like spectrum 
between 2 keV and 10 keV and a standard (not multilayer) telescope, the ~70% 
of counts are in the 2-3 keV energy range. For this reason, in the following we 
focus our attention on the performance at low energy. 

Data analysis 

The image of the photoelectron track is processed with an algorithm which 
reconstructs the absorption point as the center of gravity of the charge distri- 
bution (barycent er) and the dire c tion o f emission as that which maximizes the 
second moment ( Bellazzini et al. , 20031) . The latter is basically the direction of 



elongation of the track. However, for energies higher than ^3 keV, the initial 
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Figure 8: Examples of tracks at 2.0 keV (a) and 5.2 keV (b). The reconstructed direction of 
emission obtained by analyzing the entire track and resolving the first part are the black solid 
line and the red dashed one, respectively. The barycenter and the impact point are the black 
and the red crosses. 



part of the track with lower ionization density (energy losses are inversely pro- 
portional to the energy) and the point where photoelectron is stopped (Bragg 
peak) are resolved. The algorithm selects the pixels at the edge of the track 
with the lower signal density and calculates the impact point and the maximum 
second moment, i.e. the direction of emission of the photoelectron, considering 
only these pixels. The initial direction of photoelectrons is used to create the 
histogram of emission angles, that is the modulation curve, which shows a cos 2 
modulation if the absorbed photons are polarized while it is flat for unpolarized 
radiation. Depending on the polarization degree to be detected, the modulation 
curve must be constructed with a large number of photons, at least several tens 
of thousands, to achieve a result which is statistically significant. 

In Figure [8] we report two examples of tracks at 2.0 keV and 5.2 keV to show 
the effectiveness of our reconstruction algorithm. At low energy the track is so 
short that any substructure is blurred by diffusion during the drift in the gas 
cell. However a fraction of the original polarimetric information is still present 
in the elongation of the charge distribution, which is correlated with the initial 
direction of photoelectron emission as demonstrated by the measurement of a 
modulation factor value different from zero (see below). Instead, when tracks 
are more energetic and hence longer, the final part is clearly distinguished by the 
denser charge density. The direction reconstructed at the second step, reported 
as the red dashed line in Figure |8(b)| is evidently a better approximation of the 
actual direction of emission than that obtained at the first step as the elongation 
of the charge distribution (black solid line in Figure). As a matter of fact, the 
improvement of the modulation factor, that is the response to polarization, 
achieved passing from the first to second step is 28% at 5.2 keV and 50% at 
7.8 keV. 
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Figure 9: Systematic effects which emerge by selecting only tracks which hit less than 27 pixels 
at 2.0 keV. The 60° periodicity reflects the hexagonal pattern of the GEM, whose pitch is 
larger than usual. 

Beyond selecting spatially the incident beam and excluding the tracks with 
more than one cluster, as explained in Section G2 we remove from the analysis 
on modulation factor those events which (i) hit less than 27 pixels and (ii) have 
an asymmetry less than a threshold that is different for each energy. The first 
cut is directly related to the GEM used, which has a pitch larger than usual 
(80 /im). When primary electrons are amplified by the GEM, the charges are 
constrained in the holes and this actually samples the event on the microscopic 
structure of the GEM. If the track is so short to involve only a few holes, the 
amplified charge distribution reflects significantly the GEM pattern. This is 
hexagonal and then produces on the elongation of the track, ultimately related 
to the reconstructed direction of emission, a systematic effect with periodicity 
60°, namely the amplification makes most probable the directions aligned with 
GEM holes, which correspond to integer multiples of 60°. Pixels on the ASIC 
are arranged more finely than holes in the GEM (pitch 50 /jm vs 80 /itm) and 
removing those events which hit less than 27 pixels implies that we exclude from 
analysis tracks passed through only 10 GEM holes or less. 

The modulation curve for events which hit less than 27 pixels is reported in 
Figure [H The sharp peaks clearly mirror the hexagonal pattern of the GEM 
and significantly affect the total modulation curve because the fraction of events 
with less than 27 pixels is 10% at 2.0 keV (7.6% at 2.6 keV). 

The asymmetry of a track (or eccentricity e) is defined as the ratio between 
the maximum and minimum values of the second moment of the charge distri- 
butions M 11 , e = M^ ax / Mmin 11 . A higher value of e means that the track is 
more developed in a certain direction and then the information on the initial 
direction of emission is more preserved. On the contrary, almost round tracks 
(e ss 1) have no preferred directions and then the large part of polarimetric 
information is lost. 

We used the second cut on eccentricity to make our results comparable 
to what previously presented for the GPD filled with helium and DME by 
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Figure 10: Modulation measured for 2.0 keV polarized photons. 
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0.9 


9112 



Table 4: Modulation factor measured at different energies. Cuts are applied on the number 
of hit pixels and on the eccentricity of tracks e. The fraction of data removed (FdR), 

the reduced xi of the cos 2 fit (97 d.o.f.) and the counts after cuts are also reported. 
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Muleri et al.l ((2008). These authors reported the modulation factor with dif- 
ferent data selections aimed to optimize the quality factor /J.y/e, and, from that 
analysis, best results were obtained by removing about 25% of events. In this 
paper we follow a specular approach as for previous measurements, we change 
the cut on the asymmetry of tracks to exclude from analysis about 25% of tracks 
to maximize the quality factor. 

4-2. Comparison to Monte Carlo results and other mixtures 

An example of the modulation curve obtained at 2.0 keV, particularly sig- 
nificative because it refers to the lowest energy at which a photoelectric po- 
larimeter has ever been tested, is reported in Figure [TUJ At this energy the 
modulation factor is 13.5%, while the values measured at other energies are 
reported in Table |U 

Measured modulation factor is compared to the estimates derived from 
Monte Carlo simulations in Figure QTJ Noteworthily we achieved expected per- 
formance at low energy (even if with modest cuts), while at higher energies 
measured values are systematically lower than expected. This discrepancy is 
quite small, <10%, and actually does not impact on the sensitivity of the in- 
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strument because emission from astrophysical sources is strongly concentrated 
at low energy where expected performance is entirely confirmed. Nonetheless 
this issue is probably a signature that Monte Carlo can still be refined and the 
discrepancy will be investigated to exploit in the best way the whole energy 
range of the instrument. The measurement of possible variations of polariza- 
tion signature with energy is at the basis of the success of a X-ray polarimetry 
mission and the exploitation of the largest energy range possible, even better if 
it extends the "classical" 2-10 keV energy interval, is mandatory. 

Finally it is useful to compare the modul ation factor of the G PP presented 
in this paper with that already reported by iMuleri et all (|2008l ) for a mixture 
composed of 20% helium and 80% DME at 1 atm. Despite the fact that the 
partial pressure of DME is (almost) the same for the two gases, there are some 
differences. While at lower energies the behavior is quite similar, the DME 
mixture apparently is less effective at high energy and at 5.2 keV the quality 
factor is ~10% lower. Hence the presence of a small fraction of helium seems to 
favor reconstruction of longer tracks. 



5. Conclusion 

The characterization of the GPD presented in this paper represents one step 
forward to the optimization of the instrument as a photoelectric polarimeter. In 
place of the standard mixture composed of a small fraction of helium (20%) and 
DME (80%) at 1 atm, we tested the GPD filled with pure DME at 0.79 atm. 
Measured modulation factor successfully confirms Monte Carlo simulations at 
lower energy, while above 3.7 keV there is a discrepancy <10% with respect 
both the expected value and the sensitivity measured for the (similar) standard 
mixture. Although this issue does not significantly affect the performance of 
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DME, this work He-DME 



Energy 
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2.0 


0.135 ± 0.005 


0.058 ± 0.002 






2.6 


0.267 ± 0.004 


0.103 ± 0.001 


0.276 ± 0.014 


0.101 ± 0.005 


3.7 


0.402 ± 0.003 


0.105 ± 0.001 


0.431 ± 0.012 


0.107 ± 0.003 


4.0 


0.426 ± 0.006 


0.101 ± 0.001 






5.2 


0.496 ± 0.005 


0.082 ± 0.001 


0.545 ± 0.010 


0.091 ± 0.002 


7.8 


0.589 ± 0.010 


0.053 ± 0.001 







Table 5: Comparison of the modulation factor /t and quality factor [i\f k measured for DM E 
at 0.79 atm and a mixture He 20% and DME 80% at 1 atm presented bv lMuleri et al.l (|2008l V 

the GPD, which are mostly determined at low energy, it deserves further inves- 
tigations to exploit in the best way the whole energy range of the instrument. 
Conversely the modulation factor measured at 2.0 keV, the lowest energy ever 
presented for a photoelectric polarimeter, is relatively large, 13.5%, and encour- 
ages attempts to reduce the threshold of the instrument below this energy value. 
This could allow to address interesting scientific objectives, like the study of the 
(possibly) highly polarized thermal emission from the surface of cooling neutron 
stars. 

We also discussed, for the first time in a systematic way, the spectral capa- 
bilities of the GPD. The spectrum of absorbed photons is obtained by summing 
the charge content of hit pixels and, in principle, it could be heavily affected by 
nommiformitics in the response of the pixels. As a matter of fact, the energy 
resolution is quite good, 24% at 5.9 keV, that is very close to that of standard 
proportional counters and to the requirement (20%) defined for future space 
missions that include the GPD. The dependency with energy, beyond a \j\f~E 
contribution that is naturally explained as due to Poisson fluctuations in pro- 
duction and amplification of primary charges, shows an additional term whose 
origin will be subject of further investigation. 

In conclusion, DME appears as an interesting alternative to mixtures of he- 
lium and DME. Despite the small reduction of the modulation factor at high 
energy, performance at low energy is confirmed. Possibly even more important, 
DME provides very good spectral capabilities, which are a mandatory comple- 
ment of the polarimetric sensitivity to pursue the scientific objectives of any 
X-ray polarimetry mission. 
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